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To quantify the distribution of bar shapes in spiral galaxies, we have analysed 
113 if -band and 89 5-band galaxy images from the Ohio State University Bright 
Spiral Galaxy Survey. Parameters measuring bar shape and position along the Hubble 
sequence were obtained in each waveband. Evidence was found for a bimodality in 
the distribution of bar shape, implying that barred and unbarred galaxies are not 
just the extrema of a single distribution, and that any evolution between these two 
states must occur on a rapid timescale. Objective bar shapes measured in the H- 
band were found to be more closely related to visual classifications than 5-band 
bar strengths, as the 5-band images are somewhat compromised by localised star 
formation, especially in later-type systems. Galaxies were found to be more centrally 
concentrated in the infrared. Later type galaxies showed greater asymmetry in the 
optical than the infrared, presumably again due to localised star formation, but on 
average the bar strengths in the two bands were found to be the same. 
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1 INTRODUCTION 

For over 75 years, bright galaxies have been classified using 
Hubble's (1926) Tuning Fork. Modifications and refinements 
have been made, but the original scheme, denning "early" 
to "late" types by the size of the bulge, and the smooth- 
ness and pitch angle of the spiral arms, has proved to be a 
remarkably robust starting point. Hubble's second param- 
eter, whether or not a galaxy contains a central bar, has 
also proved a useful defining feature in describing galaxies' 
morphologies. However, the existence of galaxies with dif- 
fering degrees of "barriness" led de Vaucouleurs (1959) to 
introduce an intermediate category, SAB, between strongly- 
barred SB galaxies and unbarred SA galaxies. This finer 
gradation raises the natural question of whether galaxies re- 
ally are naturally divided into barred and unbarred systems, 
or whether some continuum exists. This issue is clearly of 
significance in trying to understand the still-disputed ori- 
gins of the bars found in a large fraction of galaxies [see, for 
example, Sellwood (1999)]. 



In an attempt to address this question, Abraham & 
Merrifield (2000) made a quantitative assessment of the mor- 
phologies of galaxies in the Frei et al. (1996) catalogue. They 
discovered that if these galaxies are plotted in a two param- 
eter "Hubble Space," with a measure of position along the 
Hubble sequence on the x axis and barriness on the y axis, 
then the galaxies naturally split into the classical bifurcated 
tuning fork seen in Hubble's original work, with the barred 
galaxies forming a remarkably tight sequence in this space. 
The bimodal nature of the barriness parameter then implies 
that barred and unbarred galaxies do form distinct popula- 
tions. However, there were two significant shortcomings in 
this analysis. 

First, the Frei et al. (1996) catalogue was not selected in 
a manner well-suited to this type of analysis. Indeed, since 
it was chosen to contain a broad range of galaxy types, one 
might well imagine that the most strongly barred galaxies 
were included to represent this type of system, while galaxies 
with no central distortion in their isophotes were selected to 
represent unbarred systems. Thus, one could imagine that 
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the criteria upon which the catalogue was selected could 
introduce exactly the type of bimodal distribution observed. 

Second, the observations for the Frei et al. (1996) cat- 
alogue were all at optical wavelengths. Light in this part of 
the spectrum can be dominated by relatively small amounts 
of recent star formation, so the perceived morphology may 
not be representative of a galaxy's underlying structure. 
Hackwell & Schweizer (1983), and more recently Block et 
al. (1994), have shown that a galaxy's optical appearance 
can be totally different from its infrared morphology, which 
is more representative of the bulk of the galaxy's stellar dis- 
tribution. Thus, even such basic properties as whether or 
not a galaxy contains a bar may be difficult to establish 
at optical wavelengths. Indeed, it is notable in Abraham & 
Merrifield's (2000) analysis that unbarred late-type galaxies, 
where contamination by recent star formation is more of an 
issue, have tended to be classified as intermediate "SAB" 
galaxies, perhaps reflecting the difficulty in classifying the 
barriness of these systems. A further indication of the poten- 
tially misleading nature of optical classifications of barriness 
comes from the simple statistic that at visible wavelengths 
only ~ 50 — 60% of galaxies are classified as barred (Sell- 
wood & Wilkinson 1993), whereas at infrared wavelengths 
the fraction is ~ 70% (Mulchaey & Regan 1997; Knapen et 
al. 2000; Eskridge et al. 2000). 

In this paper, we address these issues by carrying out 
a morphological analysis of data from the Ohio State Uni- 
versity Bright Spiral Galaxy Survey^. The galaxies in this 
sample have well defined selection criteria, and should prove 
representative of bright spiral galaxies in the local Universe, 
removing the first concern expressed above. Further, galaxy 
images are available at both near infrared and optical wave- 
lengths, so the morphological properties can be determined 
using the more robust infrared data, and tested against the 
morphologies derived in the optical. 

The remainder of the paper is laid out as follows. In Sec- 
tion H the data set is briefly described. Section [5] describes 
the measures adopted to quantify each galaxy's position in 
Hubble Space, while Section ^| presents the results of apply- 
ing this analysis. Conclusions are discussed in Section |d| 



2 THE OSU BRIGHT SPIRAL GALAXY 
SURVEY 

The sample of galaxies used in this investigation were taken 
from the Ohio State University Bright Spiral Galaxy Survey 
(OSUBSGS). 205 galaxies were observed in B,V,R,J,H, and 
K bands at five different observatories, using 1.2m - 2.4m 
telescopes. A full description of the survey can be found 
in Eskridge et al. (2002). The spiral galaxies were selected 
for the survey from Third Reference Catalogue of Bright 
Galaxies (RC3) (de Vaucouleurs et al. 1991) subject to the 
constraints that they have a B-band magnitude B < 12 and 
a diameter D < 6'. 5. 

In this paper a total of 196 -ff-band, and 166 B-band 
galaxy images were investigated. However before the classi- 
fication parameters could be measured, the foreground stars 

* The subset of the OSUB SGS used was the Early 



early 



ttype 



late 



Data Release, available o nline at http://www.astronomy.ohio- 
3tate.edu/~survey/EDR/ 



Figure 1. Distribution of galaxy Hubble types in the H and B- 
band subsamples of OSUBSGS, in comparison with the RC3. Ni 
is the normalised number of galaxies in the RC3, while rii is the 
normalised number of galaxies in the OSUBSGS subsamples, so 
for an unbiased sample rii/Ni = 1. 



in each image had to be removed. Stars were replaced with 
an average of the surrounding pixels. All images were closely 
examined for suitability and were rejected if they had too 
many foreground stars; had bright stars near the centre; or 
had large bright stars elsewhere causing saturation over a 
significant portion of the galaxy image. 

As it is almost impossible to determined whether a 
galaxy that lies close to edge-on is barred on the basis of 
photometry alone, we follow Abraham et al. (1999) and only 
use galaxies whose axis ratios imply an inclination i < 60°. 
In the B-band images, galaxies with strong dust or star for- 
mation features (e.g. NGC 2442), could not be satisfactorily 
fitted by any kind of ellipse-fitting model, and were therefore 
also rejected. In total 89 B-band and 113 //-band images 
were found to be suitable, with a total of 72 galaxies having 
images in both bands. Fig. |l| shows that the sub-sample is 
slightly deficient of very late-type systems (due to the re- 
jection of heavily dust-obscured galaxies), but is otherwise 
representative of the nearby bright spiral galaxy population 
found in the RC3. 



3 CLASSIFICATION PARAMETERS 

In order to classify the galaxies, we adopt the framework 
described in detail in Abraham & Merrifield (2000). This 
methodology defines a two-dimensional Hubble Space where 
the ^-coordinate measures position along the early-to-late 
sequence, while the y-coordinate measures, in a quantitative 
way, the degree to which a galaxy is barred. This parameter 
space provides a quantitative framework for investigating 
Hubble's Tuning Fork: if the original scheme were strictly 
valid, we might expect to find galaxies separated into one- 
dimensional barred and unbarred sequences, with weakly 
barred galaxies perhaps lying at the lower edge of the barred 
distribution. 

We assume that the distribution of light corresponds to 
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a thin disk that is intrinsically axisymmetric at large radii, 
and define the centre of galaxy as the pixel with the max- 
imum flux. A quantitative measure of asymmetry, A, was 
obtained by subtracting the image of the galaxy from a ver- 
sion that had been rotated through 180 degrees (Abraham 
et al. 1996). We define the early-to-late (x) axis using the 
central concentration parameter, C, defined in Abraham et 
al. (1994) and closely related to the parameter defined by 
Doi, Fukugita, & Okamura (1993). This parameter closely 
tracks bulge-to-disk ratios and has been shown to provide 
a quantitative substitute for more orthodox visual classifi- 
cations of position along the Hubble Sequence (Abraham et 
al. 1996). The concentration parameter is simply the flux 
ratio between an inner and a outer ellipse. The outer ellipse 
is selected using second order moments obtained by a 2 a 
cut above the sky noise, a. The inner ellipse is calculated by 
a similar moment analysis, but at a radius of only 30% of 
the outer ellipse. 

Assuming that the galaxy is intrinsically axisymmetric 
at large radii, the axis ratio of the outer ellipse can be used 
to define the galaxy inclination, which is also important for 
determining the bar shape (y) axis. Following Abraham & 
Merrifield (2000) we define the bar strength as; 



H-Band NGC3938 



B-Band NGC3938 



fbar 



arctan 
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where (^)bar is the intrinsic axis ratio of the putative bar, 
calculated from its apparent axis ratio and the galaxy's incli- 
nation using equation (2) of Abraham et al. (1999). To make 
sure that we do not miss any small high surface brightness 
or large low surface brightness bars, the whole range of radii 
is searched starting with the outer isophote (defined as for 
C) and working inwards. The inner ellipse is then defined as 
the isophote with the minimum {^)bar (he. the maximum 
value of /bar)- Figure [] shows several examples of the ellipse 
fits. 

In order to discriminate between barred and unbarred 
galaxies, Abraham et al. (1999) proposed that systems with 
(a)'bar < 0-5 could be classified as barred. This corresponds 
to fbar < 0.11, which we will use to determine the fraction 
of barred galaxies in the OSUBSGS. It should be stressed 
that this cutoff does not represent any underlying physical 
process within spiral galaxies, but merely provides a useful 
visual criterion for defining a galaxy as barred. 



4 RESULTS 

Figure |§] displays the distribution of the optical and near 
infrared samples in the two dimensional Hubble Space. A 
slight anticorrelation is observed between C and /b ar in both 
infrared and optical bands. As discussed in Abraham & Mer- 
rifield (2000), such a trend occurs naturally as the presence 
of a strong bulge (and hence large value of C) will tend to 
wash out any bar signature, whereas in galaxies with weak 
bulges even quite modest bars will still be strongly detected. 

In terms of their qualitative bar classifications, the 
galaxies are reasonably well divided in both optical and near 
infrared Hubble Space, with SA galaxies distributed across 
the bottom, SB galaxies across the top, and SAB galax- 
ies bridging the gap between the two. The bimodal split 
between barred and unbarred galaxies is far less apparent 





B-Band NGC4136 




H-Band NGC3513 



B-Band NGC3513 





Figure 2. Ellipse fits for H and B band images. The solid black 
line is the outer ellipse used to determine the galaxy inclination, 
and for calculations of C and ft ar . The white ellipse is the inner 
ellipse used for measuring f oar . The dashed lines are isophotes 
fitted to different intensity slices separated by lcr in the sky noise. 



in these figures than was found by Abraham & Merrifield 
(2000). However, there is still evidence for such a split even 
in this more objectively defined sample. The histograms 
of /bar in Fig. H show a dip in the number of galaxies at 
/bar ~ 0.2, the same value as the gap found in Abraham & 
Merrifield's (2000) analysis of the Frei sample. The dip is at a 
greater significance level (~ 95%) in the infrared data than 
in the optical, which suggests that the optical bar shapes 
are more compromised by the effects of star formation, dust 
obscuration, etc. 

This impression is borne out by the fact that there is a 
tighter grouping of the visual classifications, particularly the 
SB galaxies, in the infrared data than in the optical (despite 
the fact that the visual classifications were based on optical 
images). The 5% of cases where the visual and automated 
classifications disagree strongly were inspected individually. 
In almost all cases the discrepancy arises because the auto- 
mated fit is being compromised by bright star-forming fea- 
tures. In some cases, these features are localised regions, 
but in others (such as NGC4136, see Fig. |^) the fit is drawn 
to tightly wrapped spiral arms, which are identified as the 
strongest non-circular features, rather than the bar that we 
seek. It should be noted that for galaxies with two or more 
bars, the automated process would pick out the most ellip- 
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Figure 3. Distribution of optical and near infrared samples in "Hubble Space". SA galaxies are represented by circles, SAB galaxies by 
crosses, SB galaxies by filled circles, and unclassified galaxies by asterisks. 



Bband 



a 
CD 



0.2 0.3 



a 
CD 



Hband 



Figure 4. Histogram of bar shape distribution in (a) Optical and 
(b.) Near Infrared. 



tical bar, not necessarily the primary outer bar (although 
no such cases were found in this sample). We could have 
intervened in the bar analysis to prevent such misfits, but 
such subjective judgements are contrary to the ethos of this 
approach, as they could easily introduce subjective biases in 
the analysis. 

There are thus at least two non-physical sources of scat- 
ter in the distribution of /bar- First, the intrinsic limita- 
tions of the automated fitting process (particularly at optical 
wavelengths) mean that there is an uncertainty introduced 
in the derived value of bar shape. Second, the correlation 
between C and /bar means that a projection of this two- 
dimensional space onto the /b ar axis will tend to smear out 
intrinsically sharp features in Hubble space when viewed in 
one dimension. Hence, the true divide between barred and 
unbarred galaxies is likely to be rather more dramatic than 
we see in Fig. ^. 

In order to compare the fraction of barred galaxies in 
the optical and infrared, only the 72 galaxies for which im- 
ages in both bands were available were used. We find that 
74% of galaxies are barred in the optical, with this fraction 
increasing to 79% in the near infrared. Visual classification 
of the entire OSUBSGS sample, preformed by Eskridge et 
al. (2000) found only 72% of galaxies to be either strongly 
or weakly barred in the //-band, while only 64% were classi- 
fied as barred in the RC3 (B-band). Overall there is a good 
correlation between optical and NIR bar shape. Obscura- 
tion from dust and star-formation did cause some late-type 
galaxies to appear more barred in the NIR than the optical. 
A few galaxies did appear to be significatly more barred in 
the optical than NIR, as isophotes fitted around HII regions 
can occasionally give the appearance of a bar. 

Automated classification finds a slightly higher fraction 
of barred galaxies in both bands than visual classification. 
This could mean that smaller bars are occasionally over- 
looked when visual classifications are made. However the dif- 
ference in fractions in this case is probably due to the value 
of /bar selected to separate barred from unbarred galaxies. 
Indeed, if we increase the proposed value of fb ar < 0.11 to 
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The comparison between optical and near-infrared mor- 
phologies is also instructive. Some measures of galaxy mor- 
phology, such as central concentration and asymmetry, vary 
systematically between the optical and near infrared. How- 
ever the bar shape parameter remains approximately the 
same, albeit with greater scatter in the optical due to lo- 
calised star formation, etc. This discovery has an important 
implication for studies of bar shape as a function of redshift, 
which indicate that barred galaxies are extremely rare for 
z > 0.5 (Abraham et al. 1999; Van den Bergh et al. 2000). 
One possible explanation for this rarity is that it might be 
attributed to bandshifting effect, since the /-band images 
used in these studies correspond to B-band rest-frame emis- 
sion for z ~ 0.8. Thus, if bars were systematically weaker in 
the B band than at red wavelengths, then one would expect 
to see such a decrease in bar fraction with redshift. This 
matter is discussed in greater detail by Van den Bergh et 
al. (2002) . Since there is no such systematic variation in bar 
shape with waveband, this explanation now seems unlikely; 
it would appear that bars are intrinsically rarer at higher 
redshifts. 

These conclusions have been drawn from a relatively 
small sample of galaxies, necessarily limiting their confi- 
dence. However, with the advent of large surveys like the 
Sloan Digital Sky Survey (York et al. 2000), much larger 
samples of galaxies are becoming available. Since the mor- 
phological parameters used in this analysis are reasonably 
robust, they can be applied to rather more distant galaxies 
than traditional classifications, opening up a larger region of 
space to such studies, and hence allowing much larger sam- 
ples to be gathered. With these larger samples, objective 
morphological analysis should provide definitive answers to 
many of the questions where this study has only been able 
to scratch the surface. 
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